The interaction of dark matter cusp 
with the baryon component in disk galaxies 
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In this paper we examine the effect of the formation and evolution of the disk galaxy on the 
distribution of dark halo matter. We have made simulations of isolated dark matter (DM) halo 
and two component (DM + baryons). N-body technique was used for stellar and DM particles and 
TVD MUSCL scheme for gas-dynamic simulations. The simulations include the processes of star 
formation, stellar feedback, heating and cooling of the interstellar medium. The results of numerical 
experiments with high spatial resolution let us to conclude in two main findings. First, accounting 
of star formation and supernova feedback resolves the so-called problem of cusp in distribution of 
dark matter predicted by cosmological simulations. Second, the interaction of dark matter with 
dynamic substructures of stellar and gaseous galactic disk (e.g., spiral waves, bar) has an impact on 
the shape of the dark halo. In particular, the in-plane distribution of dark matter is more symmetric 
in runs, where the baryonic component was taken into account. 



1. INTRODUCTION 



Recent progress in understanding of the formation and 
early evolution of galaxies allows to construct a fairly 
clear picture of these processes [1-3]. According to cur- 
rent cosmological ACDM-models, gravitationally bound 
objects appeared from density fluctuations, which evolved 
hierarchically, i.e. in the «bottom-up» assembly of cos- 
mic structures. Lowest mass objects were formed first, 
and the more massive structures were formed by mergings 
and accretion of these objects. These structures mostly 
constructed from the dark matter (DM), which average 
density was approximately 6—10 times greater than the 
density of baryonic component. As a result, a large num- 
ber of gravitationally bound objects were formed in a 
wide range of masses. Baryonic matter gradually settled 
in DM potential wells. First stars (population III) and 
protogalaxies were formed due accumulation of baryonic 
gas. First stellar population has changed the chemical 
composition of the primordial gas, enriching it in metals. 

Compactness of the universe at early stages of evolu- 
tion (z = 10 — 30) caused numerous mergers of DM-halos 
with each other, forming the large-scale structures. Ap- 
proximately at z = 1 — 10 galaxies were formed. Merger 
of galaxies with similar sizes and masses (large merging) 
occurred mostly until the epoch of z ~ 1 — 2. About same 
time galaxies acquired most of the angular momentum. 
Further the galaxies, if they were not in a dense clus- 
ter, could have only small merging, i.e. merging with 
satellites of masses less than 10% of the host galaxy. 
However, despite the apparent success of the construc- 
tion of adequate models of the formation and evolution 
of galaxies, many issues remain unresolved. The essence 
of these questions is the discrepancy between the theoret- 
ical conclusions of models and observable characteristics 
of galaxies. Here are some of the problems typical for 
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cosmological models of formation and evolution of disk 
galaxies. 

• The density distribution of dark matter in the cen- 
tral regions of galaxies according to the calculations 
is much more than ocute the observable one (the 
problem of the central cusp) [4] . 

• The correlation between mass and size in models 
of galaxies is less pronounced than observations 
show [5]. 

• In the models, galaxies are generally not formed 
without a central bulge [6]. 

• Galaxies in the models only partially are balanced 
by the rotation [7]. 

• «Missing galaxy problem» : DM simulations predict 
more subhaloes (low-mass galaxies in groups) than 
number of dwarf galaxies observed [8] . 

These problems are particulary linked to insufficient nu- 
merical resolution: typical spatial scale of the baryonic 
matter evolution is of many orders of magnitude smaller 
than the spatial scale of dark halo evolution 10^ — 10^ 
kpc). Even for the galactic disk 10 kpc) there is 
no comprehensive numerical methods, adequate including 
by star formation processes, chemical enrichment of the 
interstellar medium, supernova explosions, heating and 
cooling of the gas, e.t.c, because the spatial scales of these 
processes usually do not exceed few parsecs. Therefore 
all models have to work in the approximation of «subgrid 
physics ». In any case, increasing of the spatial resolution 
and accurate treatment of stellar-gaseous component in 
the models is an important direction to improve the the- 
oretical approach. 

In this paper we considered the influence of the forma- 
tion and evolution of the disk galaxies on the dark matter 
halo distribution. A particular emphasis was made on the 
central cusp problem and on the baryonic matter influ- 
ence on the shape of a dark halo, during the disk galaxy 
formation and evolution. 
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The problem of the central cusp is well known from 
cosmological simulations of the evolution of dark matter 
and the formation of baryonic structures. The problem 
is that in numerical models the volume density p in the 
center of the DM halo tends to singularity, forming the 
so-called central cusp. But in the observations such cusps 
of the density distribution are not detected. The density 
profile in the DM-halo is described by the approximation 
p oc r" where —1.5 < a < —1. For example, for the most 
commonly used profile NFW [9]: 



(r/r,)(l + (r/r,)2) ' ^' 

where r is the distance from the halo center, po a-nd are 
the parameters which describe a particular halo, a = — 1. 

Observed density distributions of stars correspond to 
the distributions of density in the DM-halo with a > — 1 
(«core» type distribution). For comparison with obser- 
vations they commonly use data on the rotation curves 
(or velocity dispersion) in the central regions of dwarf 
galaxies [10]. Dwarf galaxies are more appropriate ob- 
jects to determine the distribution of the density of dark 
halo, which can dominate even in the central regions of 
these galaxies [11]. For more massive galaxies, observa- 
tions also confirm the «core» type density distribution of 
dark halo. For example, in [12[ detailed rotation curves 
were recovered, based on Ha kinematic survey of 36 close 
spiral galaxies. The results of decomposition of rotation 
curves of these galaxies argue in favor of quasi-isothermal 
density distribution of the dark halo (i.e. a is closer to 
than to —1 as in the NFW distribution). 

Definitely many researchers tried to find the solution 
of the cusp problem. Two approaches can be considered. 
According to first approach the cusp problem is treated 
taking into account the properties of dark matter. For 
example, in [4, 13[ it was proved that the cosmological 
random motions of «heated» DM particles in collapsing 
proto-halos. Lead to the suppression of the cusp-like den- 
sity profiles inside the forming halo, promote the forma- 
tion of DM-nucleus in galaxies and allow to explain the 
difference between observed and obtained in the numer- 
ical experiment rotation curves of galaxies. Analytical 
findings made in this approach, should be confirmed by 
the numerical N-body models, this is possible with the 
improving of the spatial resolution of the central regions 
of the halo. In [14] CDM model is considered as an ensem- 
ble of elementary particles with mixture of flavours (such 
as neutralinos; particles that may gradually «evaporate» 
from the potential well). 

According to the second, more spread approach the in- 
fluence of baryons, that have been accumulated in the 
dark matter potential well on the distribution and kine- 
matics of dark matter is considered. Following ACDM 
paradigm, various authors [15, 29[ have examined the 
«post» -evolution of dark matter under the influence of 
baryons. Baryonic matter impacts on the distribution 
of DM caused via supernovae and/or dynamic frictions 
and tidal effects. Of course, manifestation of these mech- 
anisms become substantial 1 — 2 billion years after the 



galaxy formation. In paper [16] which the linear pertur- 
bation theory and N-body simulations were used in, it 
was shown that the bar in the galactic disk may trans- 
form the dark matter central cusp (which had formed 
earlier) into the «core» type distribution. 

There are numerous models, indicating the absence of 
central symmetry in the dark matter distribution within 
the galactic disk ]17, 18]. Particularly in work ]19] on the 
basis of cosmological N-body models of galaxy formation, 
it was shown that the dark halos have a triaxial shape. 
The typical ratio for the semi-axes according to the var- 
ious data: in the plane of the disk q = h/a ^ 0.8 — 0.95, 
in the perpendicular direction s = cja ~ 0.6 — 0.7, where 
a is the semi-major axis. The parameters of the non- 
axisymmetry depend on redshift, relative mass of halo 
and other parameters of the model. Observational mani- 
festations of a triaxial dark matter distribution are given 
in several papers: distortion of the gas layer beyond 
the optical radius of the galaxy (warping) ]?, 20], non- 
spherical shape of the isophotes of X-ray emission of hot 
gas in the halo [21], dynamics of the outer galactic po- 
lar rings ]22] , observing the kinematics of hyper velocity 
stars (HVS) ]23], dynamics of tidal streams ]24, 25]. 

According to the models, the dark halo shape is more 
symmetrical in the center than at the far periphery. Ap- 
parently, on the scale lengths of the galactic disk, dark 
and baryonic matter masses are comparable and the in- 
fluence of the baryonic component on the characteristics 
of dark matter is the most signiflcant. However, large- 
scale ACDM models do not provide an opportunity to 
follow in detail the structure and transformation of the 
dark halo on the small scale lengths (within the optical 
radius of a galaxy) under the influence of the baryonic 
component. 

In this paper we have constructed a simple model of the 
formation and evolution of isolated galaxy in the gravita- 
tional fleld of dark halo taking into account the star for- 
mation, feedback (supernova explosions and the return of 
gas by normal stars) and the radiative cooling and heat- 
ing of interstellar gas. We applied this model to study 
the possible role of baryons in the blurring of the central 
cusp of the dark halo density and transformation of the 
dark halo shape. In Section 2 the formulation of problem 
is given, an algorithm, software and input physics are de- 
scribed. In section 3 main results are presented. Section 
4, contains a critical discussion of the results. The main 
findings are given in the conclusion (Section 5). 

2. INITIAL CONDITIONS AND NUMERICAL 
METHODS 

We consider evolution of dark matter halo with and 
without baryonic component. The evolution of the 
galaxy, consisting of a dark halo and the baryonic com- 
ponent was considered at 10 • 10^" yr time scale. At the 
initial moment, the matter contained within the compu- 
tational domain is made up of dark halo (lO^^A/©) and 
gas (5 • IO^Mq). In the process of evolution, the gas can 
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Figure 1. The evolution of the baryons (dark gray dots - stellar particles, black outlines — gas) in the gravitational potential of dark 

matter (light gray dots) at the times t = 10^, t = 6 ■ 10^, t = 10 ■ 10® yr 



accrete on the galaxy from the outside and also leave the 
galaxy as a result of the evolutionary dynamic processes. 
At the initial moment gas with a temperature T = 10^ 
K was evenly distributed over the halo, rotating at a rate 
of 0.5 ■ Vc{r) {Vc{r) is the circular velocity). Typically, 
in theoretical models of formation and early evolution 
of galaxies, the potential of the dark halo is defined as 
an external factor with the given density profile of dark 
matter. In our work we started the calculation from the 
quasi-isothermal distribution of DM in order to investi- 
gate of the formation of cusp-distribution of dark matter 
density in the process of dynamic evolution of the DM- 
halo. 

Numerical treatment of gas dynamics was made 
with TVD MUSCL third-order approximation numerical 
scheme in cartesian coordinate system [26]. To ensure 
conservatism of the sheme the finite-volume approxima- 
tion of the variables was used. The minmod limiter was 
used [27] in the calculations. The algorithm and software 
package based on it, were described in [28, 29] these were 
applied for solving problems of the galactic disk dynam- 
ics [30, 31]. Gas-dynamic runs were carried on a grid with 
spatial resolution of 70 pc. The dynamic of dark matter 
particles {N = 10^) and stars {N = (0-0.25) • lO*') was 
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Figure 2. The radial distribution of surface density of the stellar 
disk (t = 10 ■ 10® yr). 

calculated using N-body technique. Total gravitational 
potential of the test star- and DM-particles and gas was 
calculated by TREEcode method [32]. 

The boundary of computational domain are open. 
This enable free exchange (accretion and/or outflow) of 
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baryons, which plays an important role in the process of 
galaxy evolution. The radiative cooling processes were 
taken into account, their rate was determined according 
to [33]. 

The process of star formation was simulated as fol- 
lows. At each step of integration cells were selected, 
criteria for selection were: the surface density of gas 
a > 10^ M0/pc^ and the gas temperature T < 50 K. 
If the dynamical time in the cell Atdyn = (37r/32G'p)^/^ 
is lower than the integration step Atint, then the test 
star particle was placed into the cell. This particle can 
be represented as a star cluster (10'^ — 10''^ massive stars). 
We calculated the mass of such star particle: 



(2) 



where S^sf is the star formation efficiency, Pgas is the gas 
density, Vceii is the volume of one cell. However, we did 
not consider the initial mass function of stars in such star 
particle (star cluster) . The values of the efficiency of star 
formation ^sf were 1% and 5%. The initial velocity of 
the test particle corresponded to the gas velocity, from 
which it was formed. Kinematics of such particles was 
calculated taking into account the gravitational potential 
of the gas, dark halo and other stellar particles. 

Star clusters (test particles) during the calculation in- 
jected mass and energy into the galaxy due to two mech- 
anisms: the explosions of supernovae and stellar winds. 
Different approaches to implementation of a supernova 
explosion in the interstellar medium are used by vari- 
ous authors: the injection of thermal energy [33-35] or 
injection of the kinetic energy of the supernova rem- 
nant ]36; 37]. In our algorithm, we applied the first 
approach. At each step of the integration (Aii„t ^ 
10'^ — 10^ yr) we tested the possibility of a supernova 
explosion in the cluster. The probability of such event is 
RAtint'Ti* , where R is the number of supernova explo- 
sions per IMq per year. For Milky Way typical value is 
R = lO""'^'^ yr~^, for starburst galaxies, R reaches 

10^^^ Mq^ yr"^- We have used a control parameter Esn 
in order to vary the contribution of supernovae in these 
range of values of R. For Esn were used values 0.1 and 
0.2. When a supernova with the energy 10^^ erg explodes 
in the star cluster, this energy is partially transferred 
to the interstellar medium. Part of this energy is con- 
verted into the kinetic energy. A second mechanism for 
the transfer of matter from test particles to gas is the 
mass loss of stars. Return of gas into the ISM depends 
on stellar mass and could be significant (for example see 
[38[). We assumed that at each step of integration of 
the system of equations describing the dynamics of star 
clusters, they lose some mass, so that each particle in 
10^° years returns to the gas half of its original mass. 



3. FORMATION AND EVOLUTION OF ISOLATED 
GALAXY IN THE FIELD OF DARK HALO 

The evolution of baryons in the gravitational potential 
of dark matter shown in fig. 1. Initially gas concen- 
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Figure 3. Relative perturbation of the surface density E{r, ip) of 
the stellar disk (t = 10 ■ 10^ yr). 



trates to the midplane, where it cools down and condi- 
tions for star formation are. Herewith a thin galactic 
disk forms rapidly (~ 10^ — 10* yr). Due the gravita- 
tional instability and the dynamic interaction with the 
DM- particles galactic disc heats up and becomes much 
thicker. In the later stages of evolution the central stellar 
bulge is formed, this indicate that the formation of the 
disc is going on «outside-in». This can caused by due to 
the gaseous streams from the external region, stimulating 
star formation in the outer parts of the disk (see below). 

Stellar disk, has a typical exponental profile of the sur- 
face density a{r) ~ exp(— r/X) (see fig. 2). Furthermore, 
the stellar component density distribution demonstrate 
the presence of structure. In fig. 3 we show the relative 
perturbation of the surface density Ti(r, if): 



T,{r, if) 



(cr(r, (fi) - {a{r, (p))^ 
(cr(r, ip))^ 



(3) 



herein (...)^ is the averaging over the azimuthal angle. 
The global two-arm spiral structure with the amplitude 
about 1 — 3% is notable. It could be formed under the 
action of gravitational instability (for example see [39, 
40[) or an inhomogeneous distribution of dark matter [30, 
41]. 

The dynamics of the gas in the disk plane is shown in 
fig. 4. In addition to the regular picture of rotation, 
gas streams from external areas and the outfiows are vis- 
ible. The velocities of streams is comparable with the 
rotational speed, but the gas densities in these areas are 
smaller. 

In fig. 5 (right panel) we show the profiles of vol- 
ume DM density along the axis of rotation of the galaxy) 
for two models of the galaxy (with baryons and without 
baryons) . In a model, which includes DM only, the forma- 
tion of a central density pike is observed. There is a dif- 
ference in the radial dependence of the density profile of 
dark matter: in the calculations, taking into account the 
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Figure 4. The velocity field of the gas component in the plane 
of the galactic disl^. White isolines show the surface density of the 
gas, gray — the surface density of the stellar disk. 



baryons, density behaves more smoothly in the central 
region (fig. 5 left panel). We do not illustrate the course 
of evolution in detail, let's give only a brief description 
of the system evolution. The baryonic gas coming into 
the galaxy is accumulated in the center. If there were no 
star formation, the concentration of baryons in the center 
will increase the compression of dark matter. However, as 
seen in fig. 5, the density distribution of DM in the pres- 
ence of the baryonic component becomes «wider». The 
reason for this is the star formation that begins in the 
dense central gas and leads to the appearance of super- 
novae, which are the powerful sources of energy and gas 
from star-forming region. As is evident from fig. 5 the 
increase in value i?sw leads to more intense ejection of 
baryons from the central region and accordingly to a flat- 
ter shape of DM distribution. Ejection of gas leads to a 
redistribution of baryons in the center, which is similar 
to the «core»-type profile. The gravitational contribu- 
tion of the baryons in the central region is comparable to 
the DM. As it'seen in fig. 5 the density of baryonic mat- 
ter in the center of the galaxy is 5 — 7 times higher than 
density of the DM. So the distribution of DM is also get- 
ting similar to the «core»-type distribution. Therefore, 
our results suggest that the «normal» evolution of the 
baryonic component of galaxies leads to the evolutionary 
smoothing of the cusp, even if it was formed earlier, and 
prevents its occurrence in the later stages of evolution. 

Now consider the problem of « large-scale » influence of 
the baryonic component on the dark halo shape. The 
distribution of dark matter within the galactic disk is 
of signiflcant interest, in terms of possible observational 
manifestations. Using our model, we examined the evolu- 
tionary changes of the shape of dark halo in the neighbor- 
hood of the galactic disk under the influence of baryons. 

In fig. 6 the spatial distributions of the parameters q 
and s (the ratios of small semi-axis of the ellipsoid to the 
semi-major axis: in the galactic plane in the perpendicu- 



lar direction, respectively) are shown. In the absence of 
baryons, the deviation of the distribution of density of the 
DM from the circular one is insignificant on the periphery 
and quite noticable in the central regions. As it is evident 
from the distribution of q in the model, which includes 
DM and baryons, the distribution of DM due to inter- 
action with baryonic matter of the galaxy in the plane 
of the disk becomes more symmetric. It seems that the 
spiral structure and/or the bar arising in galactic disk, 
interacting with a dark halo leads to the decrease of non- 
axisymmetry [15]. In the direction perpendicular to the 
disk plane, DM halo becomes more oblate, this is a con- 
sequence of the deeper baryonic gravitational potential 
well, in which the dark matter accumulates. 



4. DISCUSSION 

The initial mass of baryons in the model is 5 • IO^Mq at 
the t = 10^'^ years it increased up to lO^^M© due to the 
openness of the system changed. Despite the evolutionary 
increase of the mass, the galaxy formed i our model can 
be classified as a dwarf. As already has been noted, the 
dwarf galaxies are the most suitable objects to determine 
the distribution of the dark matter density, which these 
systems is usually dominated. Anyway DM determines a 
significant share of the potential in the central regions of 
galaxies. Three reasons for this are considered. 

First, from the shallow gravitational potential well of 
dwarf galaxy, the gas leaves the galaxy more easily than 
from the nuclei of massive galaxies. In [42] it was shown 
that there is a critical mass of gas that creates a po- 
tential well, which for single supernova explosion is deep 
enough to prevent the ejection of a substantial amount 
of gas in the intergalactic medium (in the original paper 
proto-globular clusters studied, but those are similar the 
protodwarfs galaxy). The critical mass was 
it depends on the parameters (the initial density dis- 
tribution, the details of the cooling law, etc). In our 
runs supernova explosions rate does not depend on time 
{R = const) that is not suitable for realistic galaxy evo- 
lution. In the more massive galaxies where R can be the 
two orders higher than the current rate, a significant ejec- 
tion of gas may also occur (the mode of star formation 
burst). The second reason is that the low- mass DM halo 
do not have time to accumulate a substantial amount of 
gas from their surroundings. If the hierarchical merging 
of subhalos occur, relatively low value of the densities 
ratio pb/pDM remain. Only after the reaching of a cer- 
tain size (mass) of galaxy, accretion of gas from the in- 
tergalactic medium becomes comparable to the efficiency 
accumulation of mass in the process of minihalos merg- 
ing. Third, the frequent collisions of many dwarf galaxies 
lead to the loss of gas from them, and this decreases the 
ratio pb/pDM- 

The fact that the effect of smoothing the cusp of the 
DM density expressed for our model galaxy, allows to 
assume that the effect of smoothing of the cusp under the 
influence of the evolution of the baryonic component in 
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Figure 5. The radial distribution of the average volume density of dark matter is shown on left panel: the dashed line — in the model 
without baryons, the solid lines correspond to different models with baryonic component (1 — Esn = 0.1, Esp = 0.05; 2 — Esn = 0.1, 
Egp = 0.01; 3 — Egff = 0.2, Egp = 0.05). The distribution of volume density perpendicular the plane of the galaxy is shown on left 

panel: the dashed line corresponds to the density of DM in the model without baryons, the solid line with squares — the density 
distribution of baryonic matter, the solid line — the density distribution of DM in the model, with the baryon component [Esn =0.1, 

EsF = 0.01). 
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Figure 6. Distribution of the parameters of non-axisymmetrical 
density of dark matter: open squares correspond to the model 

without baryons, dark circles — the model with baryons. On the 
left — the value of q in the plane of the disk ,on the right s — 
perpendicular to the disk plane. 



val of 10 billion years. The high spatial resolution in 
three-dimensional model has allowed to identify the role 
of small-scale processes (star formation and feedback) in 
the transformation of the dark halo shape under the influ- 
ence of baryonic matter. The mechanism of the smooth- 
ing of the density cusp in the distribution of dark matter 
in the central region of disk galaxies is proposed. The 
process of star formation, is accompanied by an outflow 
of matter from the central regions of galaxy. The cusp of 
DM density distribution is smoothed due to the gravita- 
tional influence of baryons, since the density of baryons 
in the galaxy's center is higher than the density of the 
DM. The interaction of dark matter with the small-scale 
dynamical structures of the stellar disk (spiral density 
waves, star bar) leads to a symmetrization of dark mat- 
ter halo in the galactic plane {b/a ~ 0.9 — 0.95). 



the central regions of such galaxies does occur. For more 
massive galaxies, the evolution in the central regions is 
mainly determined by baryons. 

Quite recently papers appeared, confirming our con- 
clusion on the causes of smoothing of the cusp. In works 
[43; 44] ACDM simulations of the early evolution of galax- 
ies were carried out. Authors conclude that the activity 
of supernova does flatten the initially cusp-density distri- 
bution of DM. 



5. CONCLUSION 

Evolutionary model of the formation and evolution of 
isolated disk galaxy was constructed on the time inter- 
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